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' Abstract. We compare the 3D clustering of old passively-evolving and dusty star-forming z ~ 1 EROs from the 

I K20 survey. With detailed simulations of clustering, the comoving correlation length of dusty star-forming EROs is 

constrained to be less than rg ~ 2.5 Mpc. In contrast, the old EROs are much more positively correlated, with 
I ' 5.5 ro/{h~^ Mpc) ^ 16, consistent with previous claims for z ~ 1 field early-type galaxies based on analyses 

O I of ERO angular clustering. The low level of clustering of dusty star-forming EROs does not support these to 

be major mergers building up an elliptical galaxy, or typical counterparts of SCUBA sources, but it is instead 
^ ' consistent with the weak clustering of high redshift blue galaxies and of luminous local IRAS galaxies. Current 

hierarchical merging models can explain the large ro for z ^ 1 field early-type galaxies, but fail in matching their 
high number density and overall old ages. 
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1. Introduction 

Extremely red objects {R ~ K > 5, EROs hereafter) are 
providing increasingly stringent constraints on our under- 
standing of the formation of galaxies in general, via their 
spectral evolution and clustering properties. The very red 
colors of EROs are well known to be both consistent with 
old passively evolving distant {z > 0.8) elliptical galax- 
ies (e.g. Cohen et al. 1999; Spinrad et al. 1997) or dust- 
reddened starburst galaxies (e.g. Cimatti et al. 1998; Smail 
et al. 1999). Purely passive evolution of the present day 
population of elliptical galaxies is consistent with the mea- 
sured surface density of faint EROs with K ^ 17-22, while 
current renditions of the semianalytical hierarchical merg- 
ing models fail to reproduce the surface density of EROs 
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by a large factor (Daddi et al. 2000a; Smith et al. 2001; 
Firth et al. 2001). 

Recently, we have completed a relatively large deep 
survey of very red galaxies covering 700 arcmin^ (Daddi 
et al. 2000b, DOO hereafter), concluding that EROs are 
strongly clustered in projection, by an order of magnitude 
more than all galaxies at the same limits of if < 18-19.2. 
With careful attention to the measurement uncertainty 
inherent in narrow field data, Daddi et al. (2001, DDI 
hereafter) showed that the angular clustering of EROs 
implies a spatial correlation length of = 12 ± 3 
comoving Mpc, consistent with the assumption that the 
ERO population is dominated by elliptical galaxies. This 
large clustering amplitude is not at odds with recent hi- 
erarchical merging models, which require that the most 
massive galaxies are clustered more strongly than the gen- 
eral galaxy population at high z (e.g. Mo & White 1996). 
Our results on the angular and spatial clustering of EROs 
have been substantially confirmed by the Las Campanas 
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Redshift Survey data (McCarthy et al. 2001; Firth et al. 
2001; Moustakas & Somerville 2001). 

In our recent large K20 redshift survey of a flux lim- 
ited sample of ^ 500 galaxies with K < 20 (Cimatti et 
al. 2002, C02 hereafter), we obtained redshifts for a sub- 
sample of 35 EROs. For red objects with R — K > 5 and 
K < 19.2, about 1/3 were identified as old systems (con- 
sistent with being passively evolving elliptical galaxies), 
1/3 were found to be dusty starburst galaxies and 1/3 re- 
main unidentified. While the derived fraction of early-type 
galaxies, 50 ± 20%, is consistent with previous estimates 
based on morphology (Moriondo et al. 2000; Stiavelli & 
Treu 2000), C02 showed that the dusty star- forming (SF) 
objects do contribute significantly to the ERO population 
at faint magnitudes, thus complicating the interpretation 
of both ERO surface density and clustering, as measured 
in earlier analyses. In particular, given the strong inter- 
est in the clustering amplitude of early- type galaxies, it is 
important to estimate separately the clustering properties 
of the old and of the dusty-SF EROs, hence their relative 
contribution to the clustering of the whole ERO popu- 
lation. This is attempted in this letter, where we adopt 
a cosmology with V,a = 0.7, ilm = 0.3 and Hq = lOOh 
km/s/Mpc. 

2. Clustering analysis 

2.1. The sample and the diagnostic method 

Table || shows the redshifts of the EROs identified in the 
K20 survey (C02) and classified as old passively evolving 
or dusty-SF galaxies, sorted with increasing redshift and 
divided between the two survey fields (32.2 arcmin^ from 
CDFS and 19.8 arcmin^ from 0055-27). The classification 
of EROs as old galaxies is based on the detection of the 
4000A break and Call H&K absorption with undetected 
(or very weak) [OIIJA3727 emission, while objects with 
strong [OIIJA3727 emission and an absence of a distinctive 
4000A break were assigned to the dusty-SF class (see C02 
for details). 

Despite being by far the largest sample of EROs with 
identified redshifts, standard methods for evaluating the 
full two point correlation function cannot be still applied 
because of the small number of objects. Nevertheless the 
clustering properties of the old and dusty-SF samples can 
be investigated by studying the frequency of close pairs. 
This kind of approach has been applied in regimes with 
limited amount of information, e.g. to early studies of 
QSO clustering (Shaver 1984, cfr. also Hartwick & Schade 
1990), or to analyses of the arrival directions of ultra high 
energy cosmic rays (Tinyakov & Tkachev 2001), and re- 
lates to the integral under the correlation function on 
small scales, where most of the amplitude lies. 

From Table ||, it can be noted that the sample of old 
EROs contains two pairs that, within the observational 
redshift accuracy, have the same redshift {z = 0.896 in 
the 0055-27 field and z = 1.222 in the CDFS), with an ad- 
ditional object close to the second pair aX z — 1.215. On 



Table 1. ERO redshifts in the K20 survey. All but four 
EROs have K < 19.2. The redshift measurement errors 
are preliminarily estimated to be of the order of cr 100- 
200 km/s. 
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^ Objects with 19.2 < K < 20 

the other hand, the sample of dusty-SF EROs contains no 
really close pair, the closest pair having a relatively large 
redshift separation Az = 0.015 {z = 1.094 and z = 1.109 
in the CDFS, corresponding to Av ^ 4500 km/s). The two 
old ERO pairs with the same redshift have also quite small 
angular separations ( ^ 1'), implying spatial separations 
of 0.51 and 0.82 h'^ Mpc, while the two closest dusty-SF 
pairs are separated by 24 and 40 Mpc, respectively. 
The number of independent pairs in the samples is 81 for 
the dusty-SF EROs and 49 for the old EROs, thus imme- 
diately suggesting a higher intrinsic clustering amplitude 
for the old EROs. 

To assess the significance of observed pair counts we 
first generate random samples. The selection functions are 
constructed from the observed redshift distributions of the 
two ERO populations. Simulated samples were built by 
assigning at random a redshift (rounded to Az = 0.001 to 
match the data redshift measurements) extracted from the 
appropriate selection function, with sky positions within 
boundaries matching the area of each of our fields, and 
number of objects as in the relative observations (Table 
1). The resulting probability of finding by chance > 2 
pairs of old EROs within a separation < 0.82 Mpc is 
about 5 X 10^^, a clear evidence of clustering among the 
sample of old EROs. On the other hand, the probabilities 
of finding the closest dusty-SF ERO pair at < 24 
Mpc and the two closest pairs at < 40 Mpc are both 
~ 97%, consistent with purely random chance. 

2.2. Comparison to clustered samples 

We now proceed a step further and generate simulated 
samples incorporating a known 2-point clustering ampli- 
tude, in order to derive information on the clustering of 
the two classes, and to obtain meaningful estimates of the 
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Fig. 1. For each panel we show the redshift distribution 
of 100 objects extracted from our simulations, that incor- 
porates a given correlation length rg, as indicated on the 
figure, illustrating the influence of the correlation ampli- 
tude on the overall smoothness of pencil beam redshift 
surveys (see text). 

variance inherent in the pairs statistics in the sample. We 
follow the recipe described in DDI, based on the Soneira 
& Peebles (1977; 1978) prescription, allowing us to gen- 
erate many samples with a given value of tq over very 
large volumes. We adopt the canonical parameterisation 
^(r) oc r^'' with a slope of 7 = 1.8 (justified by the ob- 
served angular slope 6 = 0.8, DOO) for the 2-point corre- 
lation function and allow the amplitude to vary. 

For these simulations one has also to account for the 
redshift space distortion, which tends to decrease the num- 
bers of small scale pairs, and for the measurement error 
in the redshift. For the pairwise peculiar velocity disper- 
sion we adopt the local value of ai2 = 360 km/s (Landy, 
Szalay & Broadhurst 1998, see also Peacock et al. 2001) 
and their functional parameterisation, which is assumed 
not to evolve significantly over the redshift range of our 
data (e.g. Kauffmann et al. 1999). For the redshift error 
a = 150 km/s is adopted (cfr. Table 1), and we note that 
its contribution is small compared to the peculiar velocity 
term. To each simulated object, an error in the redshift 
measurement and a peculiar velocity is added in quadra- 
ture, chosen randomly from the appropriate distributions, 
before rounding its redshift to Az = 0.001 to match the 
data redshift measurements. 

As expected, the close pairs statistics is strongly de- 
pendent on the correlation length. For example. Fig. ^ 
shows that in the case of strong clustering almost all the 
objects reside in spikes with 2 or more objects in each 
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Fig. 2. Top panel: the cumulative distribution of pair 
separations observed for the old EROs (heavy line with 
filled circles). The horizontal error bars show the 2cr range 
estimated from our simulations with random (solid lines) 
and clustered (dotted lines, rg = 10 Mpc) realizations. 
Bottom panel: the same but for the dusty-SF EROs. This 
comparison shows that while the error on the estimate of 
the correlation length of either sample is quite broad, it 
is clear that the dusty-SF EROs as a class are completely 
inconsistent with a correlation length of order 10 Mpc, 
estimated from projected samples of EROs (DOl, Firth et 
al. 2001). 

Az — 0.001 bin, and therefore even with our small num- 
ber of objects we would expect to find a number of very 
close pairs (as indeed we do find for the old EROs). In 
fact, with the clustered samples the probability to find 
> 2 pairs within < 0.82 h^^ Mpc increases strongly with 
tq and at the la level the observed close pairs statistics 
requires the correlation length of the present sample of old 
EROs to lie in the broad range 5.5 ro/{h~^ Mpc) ;S 16. 
On the other hand, for the dusty-SF EROs, the obser- 
vation of the two closest pairs within 40 Mpc con- 
strains rg < 2.5 Mpc at the 3a confidence level. Fig. 
H summarises concisely the comparison between the frac- 
tion of observed pairs below a given scale compared with 
the random (rg = 0) and clustered (rg = 10 h^^ Mpc) 
expectations for a range of scales. 
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2.3. Spatial and angular clustering of K < 19.2 EROs 

If we assume tq < 2.5 Mpc for the observed sample 
of dusty-SF EROs, this results in an angular clustering 
amplitude A{r) < 0.002 at if - 19. We recall that EROs 
as a whole (including both old and dusty-SF objects) have 
a factor of 10 larger angular amplitude than this (DOO). 
A solid result of this analysis is therefore that the dusty- 
SF EROs cannot be the cause of the strong angular clus- 
tering of EROs reported by DOO, in agreement with the 
considerations of DOl. It is clear from our redshift survey 
that a significant fraction of EROs are weakly clustered 
dusty-SF galaxies which therefore dilutes the true angu- 
lar clustering amplitude of the early-type galaxies popula- 
tion responsible for the majority of the clustering signal. 
A detailed estimate of the amplitude of this dilution ef- 
fect would need a more precise knowledge of the relative 
fractions of both classes and a measure of the cross corre- 
lation between the two FRO species. In fact, two dusty-SF 
EROs are in close redshift pairs with old EROs (see Table 
I), with Az < 0.002 and distances within 3.2 Mpc, 
with a probability of only 2% to happen by chance. This 
is evidence of some positive cross-correlation bewteen the 
two FRO species, an intriguing result considering the dif- 
ferent physical properties of the two populations. We defer 
a discussion of this aspect as a part of the ongoing analy- 
sis of the clustering of the whole K20 sample (Daddi et al. 
2002, in preparation). The cross correlation term will tend 
to reduce the dilution effect of the dusty-SF EROs to the 
angular clustering of all EROs. In any case, although for 
z ^ 1 early-type galaxies the spatial clustering amplitude 
of ro = 12 ± 3 h^^ Mpc (derived in DOl) is more secure, 
being based on a relatively large sample, and consistent 
with the present analysis, it is likely that such amplitude 
should be revised upward in light of the findings presented 
here. 



3. Discussion 

3.1. The clustering of dusty-SF EROs 

The clustering of dusty-SF EROs is small, and maybe con- 
sistent with the values of 1 ;S 'i'o/{h~^ Mpc) ^ 2.5 mea- 
sured for star-forming galaxies at z ~ 1 (Le Fevre et al. 
1996, Carlberg et al. 1997, Hogg et al. 2000). This would 
suggest the former to be a subclass of the latter, but with 
stronger dust extinction. Locally, dusty galaxies detected 
by IRAS are also known to have a relatively weak cluster- 
ing (e.g. Saunders et al. 1992). 

The low level of clustering seems also to be at odds 
with the idea that the dusty-SF EROs are in a starburst 
phase following a major merger event, eventually expected 
to produce an elliptical galaxy, as in this case one would 
expect to find a correlation length somewhat lower than, 
but similar, to that of the ellipticals at the same redshift. 

SCUBA sub-mm selected sources are also thought to 
be dusty objects at high redshift detected by virtue of the 
emission from dust warmed by star-formation or AGN ac- 
tivity. This population is expected (Magliocchctti et al. 
2001) and tentatively observed (Scott et al. 2001) to show 
strong angular clustering at the level of ^(1°) ^ 0.01 (see 
also Ivison et al. 2000 and Chapman et al. 2001). Therefore 
our result suggest the dusty-SF EROs are a different pop- 
ulation with respect to SCUBA galaxies, with small over- 
lap, in agreement with the latter being typically fainter 
and more distant {K > 20, median redshift z ~ 2.5-3, 
Small et al. 2000, ; see also C02, Mohan et al. 2001, and 
Dannerbauer et al. 2002 for MAMBO sources). 

Finally, we note that the class of dusty-SF EROs 
could be internally inhomogeneous: some of them maybe 
spirals with moderate extinction (e.g. van Dokkum &; 
Stanford 2001) and there could also be a mixture of dust- 
enshrouded AGNs and starburst galaxies (C02). 



2.4. Analysis of systematic effects 

We tested the stability of these results with respect to 
the statistical uncertainty in the shapes of the selection 
functions, which mostly influences the numbers of widely 
separated pairs. A change in the pairwise peculiar veloc- 
ity dispersion by 20%, would result in a change of 
only about 10% for the estimated tq values, influencing 
of course the analysis of both ERO species in the same 
direction and thus leaving the result unchanged. Fig. ^ 
shows that the two closest pairs for the old EROs in our 
survey are expected at ;$ 5 h^^ Mpc separation at the 2(t 
level for rg ^ 10 h^^ Mpc, thus demonstrating that our 
result would hold correctly even if, because of redshift er- 
rors and roundings, the two closest pairs had been found 
at Az ^ 0.001-0.002. The effect of redshift errors is in fact 
negligible for the dusty-SF ERO pairs, being all of them 
at > 20 Mpc separation. Finally, we tested that the 
result is stable to variations of the color threshold at least 
up to > 4.5. 



3.2. Field z ~ 1 early-type galaxies 

Observations of samples of faint ERO galaxies have led to 
three key conclusions regarding bright early-type galax- 
ies with L ^ L*, at z ~ 1. Firstly, their space density is 
consistent with that of local luminous early-type galaxies, 
when account is made of minimal pure luminosity evo- 
lution (PLE) (C02). Secondly, spectroscopy implies age 
^ 3 Gyr for their stellar populations (assuming solar met- 
alicity, C02); and thirdly, a comoving correlation length 
To ^ 12 Mpc (this paper and DOl) has been measured 
comparable with the local value for luminous early-type 
galaxies. 

A large correlation length, tq ^ 10 Mpc, is antici- 
pated theoretically for the hierarchical merging paradigm 
for which a rapidly increasing bias is predicted for massive 
galaxies by z ~ 1 (e.g. Mo & White 1996, Moscardini et al. 
1998). Such a large correlation length is not expected for 
the PLE (galaxy conservation) scenario (DOl). However, 
also current semi-analytical renditions of the hierarchical 
models seem to be at odds with the observed results. For 
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example, the Cole et al. (2000) model predicts a comoving 
density (Fig. 1 of Benson et al. 2001) of all the z ^ 1 galax- 
ies with IO^^Mq (consistent with om- K < 19.2 selection) 
which is a full order of magnitude below the density of just 
the old EROs observed by C02. Similarly, the Kauffmann 
et al. (1999) model Q predicts a comoving density of z ^ 1 
EROs [R - K > 5, K < 19.2) that is 3(6) times lower 
than observed by C02 for old(all) EROs. In addition, in 
these models z ~ 1 galaxies qualified as field early-types 
appear to have experienced recent star-formation, while 
the present sample of old EROs is dominated by an old 
stellar population. We conclude that to our knowledge no 
semianalytical rendition of the hierarchical merging mod- 
els can yet account for all the 3 key observed properties 
of z ~ 1 field early type galaxies described above. 
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